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We consider the extent to which future imaging surveys of galaxies can distinguish between dark 
energy and modified gravity models for the origin of the cosmic acceleration. Dynamical dark energy 
models may have similar expansion rates as models of modified gravity, yet predict different growth 
of structure histories. We parameterize the cosmic expansion by the two parameters, wq and w a , 
and the linear growth rate of density fluctuations by Linder's 7, independently. Dark energy models 
generically predict 7 pb 0.55, while the DGP model 7 ~ 0.68. To determine if future imaging surveys 
can constrain 7 within 20 percent (or A7 < 0.1), we perform the Fisher matrix analysis for a weak 
lcnsing survey such as the on-going Hyper Suprime-Cam (HSC) project. Under the condition that 
the total observation time is fixed, we compute the Figure of Merit (FoM) as a function of the 
exposure time r, exp . We find that the tomography technique effectively improves the FoM, which 
has a broad peak around t cxp ~ several ~ 10 minutes; a shallow and wide survey is preferred to 
constrain the 7 parameter. While A7 < 0.1 cannot be achieved by the HSC weak-lensing survey 
alone, one can improve the constraints by combining with a follow-up spectroscopic survey like 
WFMOS and/or future CMB observations. 

PACS numbers: 



I. INTRODUCTION 

The existence of the mysterious cosmic acceleration is usually ascribed to the presence of an extra component of 
the universe with a negative pressure, known as dark energy. However, modification of the law of gravity remains as 
another interesting and equally valid possibility. One of the most elaborated examples is the DGP cosmological model 
that incorporates the self-acceleration mechanism [3, 0] without dark energy. A fundamental question in this context 
is whether it is possible to distinguish between the modified gravity and dark energy models that have an (almost) 
identical cosmic expansion history 13, Q ■ The answer to the question is inevitably dependent on the specific model 
of dark energy or modified gravity [5j. Thus we focus on the DGP model, and consider if it has any observational 
signature that can be distinguished from dark energy models with future galactic surveys. While it is pointed out 
that the DGP model has some theoretical inconsistency at a fundamental leveljl, it is still useful as an empirical 
prototype of modified gravity models, and its observational consequences are discussed [!, [9|, [13, [H| • 

The important key is the growth rate of cosmological density perturbations, which should be different in the two 
models even if they have an identical cosmic expansion history. The weak lensing power spectrum can be sensitive 
to the growth rate, while the uncertainty of the clustering bias will be the bottleneck that makes the galaxy power 
spectrum insensitive to the growth rate. 

Currently several imaging and spectroscopic surveys of galaxies are planned to unveil the origin of cosmic ac- 
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celeration via weak lensing and baryon acoustic oscillation methods. The Hyper Suprime-Cam (HSC) project is 
a fully-funded imaging survey at the Subaru telescope, which is expected to commission in 2011. An associated 
spectroscopic survey possibility, Wide-field Fiber-fed Multi-Ob ject Spectrograph (WFMOS) project, is under serious 
discussion between Subaru and Gemini observatories (see e.g. [12|, [l3[ an d references therein for other projects). 

In the present paper, we consider the extent to which future imaging and spectroscopic surveys of galaxies can 
distinguish between the DGP and dark energy models. More specifically, we empirically characterize the growth rate 
of density fluctuations adopting Linder's 7 parameter. By optimizing imaging surveys and the combination with 
redshift survey following the previous literature [lj, KLaj . we consider how we can constrain the value of 7 from HSC 
weak lensing survey and/or WFMOS baryon acoustic oscillation (BAO) survey. 

The present paper is organized as follows: In section 2, we explain our theoretical modeling: the parameterization 
of the background expansion and the modified gravity, the Fisher matrix analysis of the weak lensing power spectrum, 
and the modeling of the galaxy sample. A demonstration with the DGP model and dark energy model is also presented. 
In section 3, our result of the Fisher matrix analysis is presented. Section 4 is devoted to summary and conclusions. 
Throughout the paper, we use the units in which the speed of light is unity. 



II. THEORETICAL MODELING 



In this analysis we consider a spatially-flat universe for simplicity, consisting of baryons, cold dark matter, and 
dark energy. We ignore the dark energy clustering, and assume that the spatial fluctuations entirely originate from 
the matter component (i.e., baryons and dark matter). We further model that the cosmic expansion history effectively 
follows the universe with the matter density parameter O m and the dark energy parameter 1 — f2 m : 
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(2.1) 



where Hq — lOO/i km s 1 Mpc 1 is the Hubble constant, a is the cosmic scale factor, and wq and w a are constants 
parameterizing the equation of state of dark energy (ltH [l7l [l8j : 



p/p = w(a) = Wq + w a (l - a). 



(2.2) 



Note that we use equation (|2.ip even in the DGP model that does not have dark energy at all by approximating 
its cosmic expansion law with the two parameters wo and w a - In this case, they do not have any relations to dark 
energy in reality, but it is already shown that such an empirical description provides a reasonable approximation to 
the cosmic expansion in the DGP model. For definiteness, the expansion in the DGP has the effective equation of 
state (e.g., pj|) 



where 
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(2.4) 



The cosmic expansion in the DGP model is well approximated by the dark energy model with effective equation of 
state with wq = —0.78 and w a = 0.32 as long as f2 ro ~ 0.27. The parameterization gives the distance redshift relation 
within 0.5 % out to the redshift 2 [H. 



A. Linder's 7 parameter 



According to refs. [19|, [20|, |21| , the linear growth factor in the DGP and dark energy models is well approximately 
expressed by 



oc exp 



a da' 
21 (fl^o') 7 - 1) 



(2.5) 



In this description, the constant parameter 7 characterizes the gravity force model, i.e., the Poisson equation. 
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The dark energy models with the effective equation of state f|2.2[) within the general relativity are well approximated 

by 

7 = 0.55 + 0.05[l + w(z = l)] (w>-l), (2.6) 
7 = 0.55 + 0.02[l + w(z = l)] (w < -1). (2.7) 

This formula reproduce the exact linear growth factor within 0.3% (0.5%) for —1.2 < w < —0.8 (—1.5 < w < —0.5). 
Therefore 7 in dark energy models takes the value 7 = 0.54 — 0.56 for —1.2 < w < —0.8 (l9l. I20I w\ . 

On the other hand, in the DGP model, the Poisson equation is modified in the linear regime. Then 7 takes a 
different value from that of the dark energy model even if the background expansion is same (i.e. if wo and w a are 



same). Ref. [2l| found that in the DGP model 7 = 0.68 is an excellent approximation for the evolution of the growth 
factor and that 7 varies by only 2 % into the past. 

The point here is that a dark energy model mimicking the cosmic expansion history of the DGP model predicts a 
different linear growth rate by A7 ~ 0.1. In what follows, therefore, we employ equations (|2.1[) to (|2.4|) to describe the 
expansion history and the growth of density fluctuations, which empirically describe both the DGP and dark energy 
models, and ask if it is possible to achieve the accuracy of A7 ^0.1 by optimizing future surveys of galaxies. 



B. Weak lensing power spectrum and Fisher matrix 

The optimization of imaging surveys is based on the weak lensing tomography method (see e.g., [22l. [23l. I24I |25| ) . 
In this methodology, one divides the entire galaxy samples in several different redshift bins according to the weight 
factor Wi(z(x)) for the i-th redshift bin: 

NiJ max{zitZ ) dz' \ x{z)J 

where dN / dz denotes the differential number count of galaxies with respect to redshift per unit solid angle (see below 
for details), x( z ) is the radial comoving distance at z, 

r — - 1 r dz> 

and 

[ z ' +1 , ,dN(z') 

Ni = i dz ^tr (2 - 10) 

is the total number of galaxies in the i-th. redshift bin. While imaging surveys provide photometric redshifts alone from 
the multi-band photometry, instead of spectroscopic redshifts, for galaxies, it is known that the lensing tomography 
works even with relatively crude redshift information. 

Assuming that the anisotropic stress is negligible, the cosmic shear power spectrum is given as: 

P(ij)(l) = / d X WMx))WMx)) ^maTs 11110 " (* - ~,*(X)) , (2-11) 

where -Pmass hncar z ) is the nonlinear mass power spectrum at the redshift z, k is the wave number of the three 
dimensional coordinates, I is the wave number of the two dimension corresponding to the angular coordinates, a is the 
scale factor normalized to unit at the redshift z = 0. We compute P^J£jJf near (k, z) adopting the Peacock and Dodds 
formula [26j]. 

The covariance matrix for P^(l) is approximately given by 

Cov[p (ii) (0,WO] - {2l + ^ /sky [i$o(Qffi)(Q + p$W)p&)(i)] 

= 5u'Cov {ijKmn) (l), (2.12) 

where we define 

Pm(l)=P(ii)(l) + ^> (2-13) 
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/sky is the fraction of the survey area, and av is the rms value of the intrinsic ellipticity of randomly oriented galaxies, 
for which we adopt a e = 0.4 (see e.g., [H, [H, [24[ ) ■ 
Finally the Fisher matrix is estimated as 

, \- \- dP fa) (Q i , n dP (mn) (l) 

F <x0-ls 1^ g 0a Cov fe)(mn)W g 9 p ■ I 2 - 14 ) 

l (ij)(mri) 

where 9 a denote a set of parameters in the theoretical modeling. To be more specific, we consider 7 parameters, 7, 
u>o 7 w<i, Q-m, (the fluctuation amplitude at 8/i _1 Mpc), h, and n s (the primordial spectral index of matter power 
spectrum), assuming the other cosmological parameters are determined from independent cosmological data analysis. 

We adopt the range of 10 < I < 10 4 x (Ng/35/rib) 1 ^ 2 for the sum of I, where N g is the number density of galaxy per 
unit solid angle (see next subsection). We define the 3 dimensional Figure of Merit by the reciprocal of the volume 
of the error ellipsoid enclosing the 1 sigma confidence limit in the {7, wo,w a } space, marginalizing the Fisher matrix 
over the other parameters. Similarly, the 2 dimensional Figure of Merit is the reciprocal of the surface of the error 
ellipse enclosing the 1 sigma confidence limit in the {wo,w a } plane with 7 fixed. 



C. Modeling galaxy sample 

We assume the following form of the redshift distribution of the galaxy sample per unit solid angle 

dN N g (3 



dz z° +1 T{{a + 1)1(3) Z 6XP 



(2.15) 



where a, (3, and zq are the parameters, and N g = J dzdN/dz. The mean redshift may be determined by 



*m = 4rf^= w, ((a + f/5 (2-16) 
N g J dz T((a + 1)/P) V ' 

We assume that N g and z m is related to the exposure time i exp as, following the reference 

/ , \ 0.067 

*» = 0.9(-f=£-) , (2.17) 
\30 mm. J 

( t \ a44 

N = 35 ( tc *P arcmin.- 2 (2.18) 

\ 30 mm. / 

The mean redshift z m changes from 0.72 to 1.1, and N g does from 7.8 to 163, as the exposure time t oxp changes from 
1 minute to 10 3 minutes. In the reference [T3 |. a = 2 and (3 = 1.5 are adopted. However, in the present paper, we 
adopt a = 0.5 and (3 — 3. 

In order to check the validity of our mock galaxy samples, we show in Figure [T] the two cases of a = 0.5 and 
(3 = 3 (dotted curve), and a = 2 and (3 = 1.5 (dashed curve), for exposure times of £ cxp = 1, 5, 10, 30, 45 minutes 
(from bottom to top respectively). The solid curves show the real redshift histograms, for the corresponding iband 
magnitude limits, taken from the CFHT photometric redshift data of [231 . These photo-z's were calibrated using the 
VVDS spectroscopy and are reliable to i ~ 25 which is sufficient for this study (see [27]). The relationship between 
magnitude limit and exposure time was scaled from the published Subaru Suprime-Cam data of [28|. These data 
are shown in Table I for the i, g, r, z passbands. Denoting the exposure time for the i band by t exp , the exposure 
time for g band is about t cxpg = 3 x t cxp . Similarly, t cxPr = 1.2 x t cxp for r band, and t cxPz = 0.3 x t CKp for z band, 
respectively. 

The total survey area can be expressed as 

Area = , f ^ld of View\ 2 2^ 

where we assume that the Field of View of 1.5 degree, the total observation time r to t a i is fixed as 800 hours, and the 
overhead time is modeled by a constant, t op = 5 minutes, plus a fraction (10%) of the exposure time Y]^ i C xpj for one 
field of view. 

We consider the cases the tomography is used, which we denote by n\, = 2, ri\, = 3 and rib = 4. Here rib denotes the 
number of the redshift bin. In the case rib = 2, the sample is divided into the two subsamples in the range 0.05 < z < z m 
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JAB limit 


i(S/N = 10) 


g(S/N = 5) 


r(S/N = 5) 


z(S/N = 5) 


22.97 




1 mins. 


3 mins. 


1.1 mins. 


0.3 mins. 


23.84 




5 mins. 


15 mins. 


7 mins. 


1.4 mins. 


24.22 


10 mins. 


30 mins. 


12 mins. 


3.5 mins. 


24.81 


30 mins. 


90 mins. 


34 mins. 


8.1 mins. 


25.04 


45 mins. 


130 mins. 


50 mins. 


13 mins. 


TABLE I: Exposure time for the bands, i,g,r,z. 


Sub — sample 


n b = 1 


n b = 2 


n b = 3 


n b = 4 


choice of band 


i 


i, r 


g,r,i,z 


g,r,i, z 




iexp 


2.2 X £exp 


5.5 X texp 


5.5 X texp 


redshift bins 




0.05 < z < 2.5 


0.05 < 2 < z m 


0.05 < z < 3« m /4 


0.05 <z<0.6xz m 








z m < z < 2.5 


32 m /4 < Z < 5z m /4: 


0.6 X Z m < Z < Z m 










5,z m /4 < z < 2.5 


Z m < Z < 1.4 X Z m 












1.4 X z m < z < 2.5 



TABLE II: Assumption on the subsample and measurement 



and z m < z < 2.5, while in the case rib = 3, we consider the three subsample 0.05 < z < 3z m /4:, 3z m /4 < z < 5z m /4 
and 5z m /4 < z < 2.5. In the case n b — 4, we consider the four subsample 0.05 < z < 0.6 x z m , 0.6 x z m < z < z m , 
z m < z < 1.4 X z m , and 1.4 x z m < z < 2.5 (see also Table II). We also consider the case the tomography is not used, 
which we denote by rif, = 1, for which we don't take into account how to obtain dN/dz, instead assuming that dN/dz 
is obtained by some method. 

We assume that the subsample of rib = 2 is constructed by the two band, r and i, observation, given that the 
strategy proposed in [29[ is successful. The cases rib = 3 and rib = 4 are constructed by the 4 band g, i, r, z, 
observation, assuming that the conventional photo-z is successful. The case rib = 1 is based on the i band observation. 
We assume that 90% galaxies of i band measurements dN/dz can be used as the subsample, in the case rib = 2, 3, 4. 

We use iexp to represent the i band exposure time for one field of view, then we assume Y) j i C xpj = 5.5 x i cxp for 
the cases rib = 3, 4, Y]^ i GX pj = 2.2 x t oxp for the case rib = 2, and J"^- t cxp ^ = < cxp for the case rib = 1, respectively. 

Figure [2] shows the resultant total survey area, and the total number of galaxies as function of the i band exposure 
time texp, for the cases, rib — 1, 2, 3 and 4. 

D. DGP model 

Here we demonstrate the weak lens power spectrum of the dark energy model and the DGP model with the same 
cosmic expansion. The linear perturbation theory in the DGP model has been extensively worked out by [ll[. While 
more recently Koyama and Silva studied nonlinear evolution of density fluctuations in the DGP model [30(, the 
nonlinear nature of the gravity in the DGP model is still an unsolved problem. Therefore we adopt an empirical 
modeling of the nonline ar g rowth combining the Peacock-Dodds nonlinear fitting formula [26[ and the linear growth 
rate in the DGP model [ll| . As a result, our predictions below may be inaccurate on nonlinear scales, but our main 
conclusions concerning the optimization strategy would be unlikely to be sensitive to this approximation. 

Figure [3] shows the weak shear power spectrum of the spatially flat DGP model and the dark energy model with 
the same background expansion. The cosmological parameters of both of the models are the same (fi m = 0.27, 
fib = 0.044, h — 0.72, erg = 0-8, and the spectral index n s — 0.95). To realize the same cosmic expansion history, the 
effective equation of state parameter of the dark energy is chosen as w(z) = —0.78 + 0.32z/(l + z), as mentioned in 
Section 2.1. A similar computation has been already considered by [3l|, but our present work differs in that we use 
the Peacock & Dodds formula and that we assume a rather shallow sample of galaxies. Because the Poisson equation 
of the DGP model is modified, then the difference comes from the growth rate. In this figure we assume 30 minutes 
exposure time of rib = 1 ■ The theoretical curves and the errors bar depend on the survey sample, but we might expect 
that the two curves could be distinguished. In the next section, we examine the capability of the differentiation. 
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III. RESULTS 



In this section, we present our optimization analyses for the HSC weak lensing survey. Specifically, we fix the total 
observation time of the HSC survey, T tot as 800 hours, and adopt a model of the background galaxy sample described 
in Sec. IIC for the HSC survey; in particular the mean redshift of galaxies z m and their surface number density JV g 
are given by eqs. (|2. 1 7|) and (|2.18[) as a function of the exposure time t exp . In this section , we also present results in 
combination with a spectroscopic survey, the WFMOS BAO survey, which will be limited by a total observation time 
(see [IH for discussion of the optimization under this condition). Note that we also assume that the WFMOS survey 
is limited by the total survey area of the HSC imaging survey. Namely, the survey area of the WFMOS survey must 
be less than or equal to that of the HSC survey, as the HSC survey is acting as a photometric source catalogue for the 
WFMOS spectroscopic survey. So, for the WFMOS survey, we fix the same survey area as the HSC imaging su rvey 



equation (|2.19|) . and the redshift range of galaxies 0.8 < z < 1.4 with the number density n — 4 x 10 4 ft, 3 Mpc 3 1 



which is a set of optimized survey parameters for the spectroscopic survey. 

Figure [4] shows the Figure of Merit (FoM) of the 3 dimension (3D) of {7, wq, w a }, as function of the exposure time, 
^exp- The 3D FoM, the reciprocal of the volume of the la error ellipsoid in the {7, Wq, w a } space, is computed by 
marginalizing the Fisher matrix of the 7 parameters over O m , ag,h, and n s , with a fixed value for the baryon density, 
n b = 0.044. 

The lensing tomography method with nj, = 2 significantly improves the 3D FoM, and continues to do so with 
increasing rib for t cxp ^ lOmins. The peak of the FoM systematically shifts to the shorter exposure time with larger 
rib, while the peak profile is fairly broad. With increasing n&, more information of redshift evolution of structure can 
be obtained. Similarly, as t cxp increases, more information of smaller structure can be obtained. However, these are 
offset by decrease in total survey area. Namely, observation of more bands and longer exposures consume observation 
time, and the total survey area becomes smaller. This decreases the FoM. 

For comparison, we plot in Figure[S]the 2D FoM, the reciprocal of the area of the la error ellipse in the {wo,w a } 
plane, evaluated by marginalizing the Fisher matrix of the 6 parameters (wo,w a ,Q m ,ag, h, and n s ) with Qb — 0.044 
and 7 = 0.55 fixed. One can find the similar features as those of the 3D FoM. This figure suggests the three redshift 
bin is enough to constrain wq and w a and that the peak of FoM is located around t cxp w 10 minutes, and the peak 
profile is very broad. The FoM of the case rib = 2 is larger than that of n b — 3,4. This indicates that observation 
of larger survey area with small number of bands (rib — 2) can be useful for the dark energy constraints, though an 
accurate-photo-z strategy is required. 

Figure [5] shows the la error on 7 as a function of t cxp , which is estimated by marginalizing the Fisher matrix of 
the 7 parameters, 7, WQ,w a , Q m , &8, h and n s , over the parameters other than 7. The curve shows the error from the 
weak lensing power spectrum adopting a proposed survey with HSC; A7 ss 0.3(1) can be achieved with (without) 
tomography. The result indicates that the weak lensing survey alone cannot reach the accuracy of A7 = 0.1 that is 
required to distinguish between the DGP and dark energy models. 

The uncertainty in 7 can be significantly (more than a factor of three) reduced by combining the baryon oscillation 
features from the WFMOS survey (Figure [7]). In modeling the galaxy power spectrum of the redshift survey, we 
simply considered the linear theory specified by the 9 parameters 7, wq, w a , f2 m , as, h, n s , bo and po, where 60 and po 
are the parameters for the bias model, for which we adopted the scale independent bias model with the form 



Here we assumed the target parameters 60 = 1-38 and po = 1. For the theoretical modeling of the galaxy power 
spectrum and the computation of the Fisher matrix, the range of the wavenumber 0.01 /iMpc _1 < k < 0.2 hMpc^ 1 is 
included, (see Appendix for details). 

From Figures [5] and the error of 7 has a minimum of t exp between several minutes and 100 minutes, depending 
on the strategy. For the weak lensing survey (HSC) alone, the tomography technique is very effective in reducing 
the error, and the result is fairly insensitive to the the choice of t exp . An additional spectroscopic survey (WFMOS) 
significantly reduces the error. In this case, shallow surveys with i exp < 10 minutes provide the minimum error for 7. 
Especially, the case rib = 1 and rib — 2 is significantly improved by the combination. This behaviour is understood 
as follows. We assume the total observation time of the WFMOS survey is not fixed, while adopting the same survey 
area as the HSC survey. Then, in these figures, the cases rib = 1 and rib =2 assumes larger survey area for the redshift 
survey than that of the cases rib = 3 and rib = 4. However, note that the minimum is located around the several 
minutes of the exposure time even for the case rib — 3 and 4. Therefore, when considering the combination with the 
redshift survey, wider and shallower surveys are indeed prefered. 

Now we are in a position to answer the question: is it possible to distinguish between the DGP and dark energy 
models ? For that purpose, A7 ;$ 0.1 is required. Figure [5] plots the 1 sigma error as a function of the total observation 
time Ttotai, where we adopt t exp — 10 minutes and n& = 4 (dash-dotted curve) and rib = 2 (dashed curve). The thin 
curve is the result of the weak lensing survey alone, while the thick curve is the result combined with the redshift 




b(z) = 1 + (b - 1)(1 + z) 



(3.1) 
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survey. Note that A7 is in proportion to T to */ 2 . Figure [8] suggests that the HSC survey alone may reach A7 < 0.1 
with Zt t = 10 4 hours, the combination with the WFMOS survey may do so with T tot = 10 3 hours if we put a prior 
constraint on f2(,- 

Finally in this section, let us consider other impact that the HSC survey may present as a test of modified gravity 
models. The dash-dotted curves in Figure [9ja) show the 1, 2 and 3-sigma confidence contours (going from the 
innermost outward) in the wq — w a plane, by marginalizing the Fisher matrix of the 7 parameters, 7, Wq, w a , £l m , 08> h 
and n s , over the parameters other than wq and w a . Here the constraint from future Planck survey is taken into 
account by including the prior constraints Afl m = 0.035, Aer 8 = 0.04, Aw = 0.32, Aw„ = 1, An s — 0.0035 [32l |. 
Here the target parameters are same as those of the ACDM model in Figure |4j and we fixed rib = 4 and t cxp = 10 
minutes. Note that the point of the DGP model (wo,w a ) = (—0.78,0.32) is marked, and is almost near the 2 sigma 
curve. This means that the HSC can distinguish between the DGP model and the ACDM model at the 2 sigma level 
by including future constraint by the observation of the cosmic microwave background anisotropy. Here, we fixed the 
total observation time as 800 hours, then the constraint can be improved when the total observation time is longer. 
The solid curve is the combination with the WFMOS survey, which also shows the significant improvement of the 
constraint. Similarly, figure^b) show thel, 2 and 3-sigma confidence contours in the wq — 7 plane, by marginalizing 
the Fisher matrix over the parameters other than wq and 7. The point of the DGP model (100,7) — (—0.78,0.68) 
is marked. With this figure, the constraint is at the 1 sigma level. Then we can not clearly distinguish between the 
DGP model and the ACDM model with this plot. These features reflect how the shear power spectrum is sensitive 
to the parameters. This suggests the choice of a projection is important for distinguishing between these models. 

IV. SUMMARY AND CONCLUSIONS 

In this paper, we investigated optimization of a weak lensing survey for the dark energy, and how such a survey 
might be used for testing modification of the theory of gravity. By introducing a simple model of the survey sample 
as a function of the exposure time for one band of one field of view, we investigated how the FoM and the constraint 
on Linder's 7 parameter depend on the exposure time and the number of passbands. To optimize the survey to probe 
probe modifications of gravity, we considered a Figure of Merit in the space {7, wq, w a } as well as in the familiar 2D 
plane {u> , w a }. We obtained the following results: 1) The peak of the FoM is located at t cxp ~ several ~ 10 minutes 
for rib = 2,3,4, though the peak profile is very broad. 2) The tomography technique improves the FoM effectively 
when including the parameter 7. 3) The combination with the redshift survey like the WFMOS BAO survey improves 
the error on the parameter 7. 4) The shallow and wide survey is advantageous for the tomography, and has potential 
when taking combination with the redshift survey into account. 5) The HSC weak lensing survey by itself is not 
sufficient for distinguishing between the DGP model and a dark energy model with the same background expansion, 
but it will be able to distinguish between the DGP and ACDM at the 2 sigma level by including the prior constraint 
from future CMB observation. 

We assumed a very simplified model of the survey galaxy sample, and the error in the photometric red- 
shift measurement is not taken into account. Also we assumed that the weak lensing power spectrum of the 
10 < I < 10 4 (Ng/35/nb) 1 / 2 can be used. Further investigation is needed including the modeling of the galaxy 
sample and the error in measuring the photometric redshift. In the present paper, we assumed the spatially flat 
universe. In general, since the lensing power spectrum is not very sensitive to the curvature of the universe, then the 
inclusion of the other parameter will degrade the constraint [33| . 
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APPENDIX A: MODELING OF THE REDSHIFT SURVEY POWER SPECTRUM 



Here we briefly review the power spectrum and the Fisher matrix formula for a galaxy redshift survey [34 , |35| , 
adopted in the present paper. Here we assume a measurement of the multipole power spectrum Vi(k) (I — 0, 2) from 
the galaxy redshift survey, which we theoretically model as 

~ m 1 f, I gMs) 2 V?(s, k, H?P(k, /x, s)Ci {jj) 

V ^=2j^ Jrfs'nW(s',fc,M) 2 ' (A1) 

where s is the coordinate of the redshift space, n(s) is the mean number density per unit volume, ip(s,k,[j,) is the 
weight factor, Ci(fj) is the Lenegdre polynomial, fi is the directional cosine between k and s, and P(fc,/i, s[z]) is the 
power spectrum at the redshift z, which is modeled as 



n/, r n s ( z ) 2 ds ( z ) o ( , ds ( z ) , n 2 S ( Z ) 



with 



Pgai.{qn,q±,z) = b{zf 



1 



(A2) 

d In Di (z) /dhia(z) q n 



b(z) q 



P^sT(q,z) (A3) 



where q 2 = q n 2 + q ± 2 , Pmass ar (9> z ) i s the linear mass power spectrum at the redshift z. The comoving distance \[ z ] 
is given by 

1 r dz ' 

X(z,n m ,w ,w a ) - — J q -j== = = _ _ = =__=__ , (A4) 



as given in equation (|2.9[) . For our fiducial model we adopt the flat ACDM model with fi m = 0.27. Thus, our fiducial 
model is s(z) — x{z, 0.27, — 1, 0). In the modeling of the bias, we consider the scale independent bias model in the 
form, Eq.flSZQ). 

The variance of Vi(k) is given by 



(2tt 



.3 



AVl{k)2 = 2 i^t S ' (S ' A:) ' (A5) 



where A 14- denotes the volume of the shell in the Fourier space, and we have defined 



1 f fd S n( S ) 4 i>( S ,k,ri 4 [P{k,»,s) + l/n( S )Y[£ l (ri? 



Qi(k)--Jdfi [/ds'n(s') 2 V(s',fc,M) 2 ] 2 ■ (A6) 



Then, we may evaluate the fisher matrix by 

1=0,2 -"Smin 

In the present paper, we adopt n(s[z]) = 4 x 10~ 4 /i 3 Mpc~ 3 and ip(s, fc,/z) = 1. 
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FIG. 1: dN/dz as function of the exposure time, a — 2, (3 = 1.5 with zo = Z m /1.41 (dashed curve), and a = 0.5, (3 = 3 with 
zq = Zm/0.64 (dotted curve), respectively, for the fitting function of the form (12. 15|) . for the exposure time t OX p = 1, 5, 10, 30, 45 
minutes from bottom to top. The solid curve shows the corresponding CFHT LS photo-z i band data. 
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FIG. 2: The total survey area (thick), and the total number of the galaxies (thin) as function of the i band exposure time t, 
for the case rib — 1, 2, 3 and 4. 
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FIG. 3: The dark (black) curve is the weak lensing power spectrum of the dark energy model with the cosmological parameter, 
fim = 0.27, Qb = 0.044, h — 0.72, erg = 0.8, n s — 0.95, and the equation of state parameter of the dark energy wo = —0.78, 
w a = 0.32, while the bright (red) curve is the flat DGP model of the same cosmological parameters. Here we assume the HSC 
like survey with £ exp = 30 minutes of the case Uf, = l(see section 2 for details). 
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FIG. 4: Three dimensional (3D) FoM in {7, u>o,u> a } as function of the i band exposure time, which is obtained from the Fisher 
matrix of the 7 parameters 7, wo, Wa, fi m , erg, h, and n 3 , Here the target parameter is 7 = 0.55, wo = — 1, Wa = 0, Q m ~ 0.27, 
cr 8 = 0.8, h = 0.72, n s = 0.95. The other parameter, Q b = 0.044 is fixed. 
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FIG. 5: Two dimensional FoM in {wo,w a } from the Fisher matrix of the 6 parameters wo, w a , fim, crs, h. Here the fiducial 
model is ACDM, with wo = — 1, w a = 0, Q m ~ 0.27, as = 0.8, h — 0.72, n s = 0.95. The other parameters, 7 = 0.55 and 
Qb = 0.044 are fixed. 
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FIG. 6: 1 sigma error in measuring 7 as function of the exposure time, obtained by marginalizing the Fisher matrix of the 
7 parameters 7, ivq, w a , fi m , eg, h, and n s , over the parameters other than 7. The target parameters is the same as those of 
Figure 0] 
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FIG. 7: Same as figure [6] but the considering the case of the weak lensing power spectrum combined with the galaxy power 
spectrum of the redshift survey. 



17 



10.00 



0.01 



1 1 1 1 1 1 1 1 1 1 


i i i i i i i | i i i i i i i i_ 

Thin HSC : 
Thick HSC + WFMOS : 


X - ' -v 

— »X , * ^ 

s ^ *••-. 




*» x, 

*X * 
** 

"X 

X. 


"X , ^ 

x» " ^ 
x» , * 


- n b -4 


•x «x , _ 

--:>:.:; 

x, 

X* 

x» 

Xff 


n b =2 

i i i i i i i i_l i 





10 100 1000 10000 

Ttotai (hours) 



FIG. 8: 1 sigma error on 7 as function of the total observation time. Here we fixed t eX p = 10 minutes and nt = 4(dash-dotted 
curve) and rib — 2 (dashed curve) . The thin curve is the result with the 7 parameters of the Fisher matrix for the lensing power 
spectrum, but the thick curve is the constraint from the combined weak lensing power spectrum and galaxy power spectrum 
(from a redshift survey). 
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FIG. 9: (a, Left) The 1, 2 and 3-sigma contours in the wo — w a plane. The dash-dotted curve is the result with the 7 
parameters of the Fisher matrix for the lensing power spectrum and the Planck prior constraint, and the solid curve is these 
constraints combined with the galaxy power spectrum from a redshift survey. The target model is the A CDM model, then 
(wq, w a ) = ( — 1, 0), and the mark (wo,w a ) = (—0.78, 0.32) is the DGP model. Here we fixed n b = 4, t cxp = 10 minutes, and the 
total observation time, 800 hours, (b, Right) Same as (a), but with the contours in the wo — 7 plane from marginalizing the 
Fisher matrix of the 7 parameters over all other parameters. The target model is the ACDM model, then (wo, 7) = (— 1, 0.55), 
and the mark (w ,7) = (-0.78, 0.68) is the DGP model. 



